Microscopic mechanism for mechanical polishing of diamond (110) surfaces 
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Mechanically induced degradation of diamond, as occurs during polishing, is studied using total- 
energy pseudopotential calculations. The strong asymmetry in the rate of polishing between different 
directions on the diamond (110) surface is explained in terms of an atomistic mechanism for nano- 
groove formation. The post-polishing surface morphology and the nature of the polishing residue 
predicted by this mechanism are consistent with experimental evidence. 



PACS numbers: 81.40.Pq, 62.20.Fe, 71.15.Nc 

The extreme mechanical properties of diamond are of 
both practical and academic interest, and have been the 
subject of much recent research [Q. Diamond, the hard- 
est known natural material, is extremely wear resistant 
and conventional polishing requires the abrasive action 
of diamond powder to wear down the surface. In spite of 
extensive experimental investigation, very little is known 
about the dominant microscopic mechanisms which occur 
during polishing, where the conditions are such that both 
mechanical and chemical processes can play a role. In 
this Letter we describe the results of ab-initio quantum - 
mechanical simulations of particular wear processes on 
the (110) surface. These simulations suggest that the 
atomistic mechanism of nano-grooving may be respon- 
sible for the main features of diamond polishing on this 
surface M. 

The most striking feature of diamond polishing on the 
(110) surface is the strong asymmetry in the rates of 
material removal between different directions. Ratios 
of greater than 10:1 between the (001) (soft) direction 
and the (110) (hard) direction are typically observed ||. 
Scanning tunneling microscopy studies p| of the post- 
polishing surface morphology also show significant differ- 
ences between the directions. In the soft direction there 
are grooves which are flat bottomed on the atomic scale 
and which are absent in the hard direction. For the soft 
direction, the absence of fractures or dislocations indi- 
cates that the regime in which polishing occurs is below 
the critical stress. Furthermore the size of clusters in 
the residue, and the small amount of material removed 
per traversal in sliding indenter experiments B , suggest a 
mechanism in which only a few carbon atoms are removed 
in each tribological event. In other respects the hard and 
soft polishing directions are more similar, for example, 
the increase in polishing rate in hydrogen and oxygen rich 
environments B and the weak temperature dependence 
of the rate ||. Studies of the polishing residue show 
both similarities and differences between the directions. 
Electron energy loss spectroscopy |j] and transmission 
electron microscopy fjj studies of the polishing residue 
indicate that it consists mainly of amorphous clusters, 



although nanometer scale regions with a structure sim- 
ilar to graphite are also observed. When polishing in 
the hard direction, additional small diamonds (~1 nm 3 ) 
are occasionally observed in the residue, although it is 
unclear whether these originated in the polishing dia- 
mond or the polished diamond, and there are also fewer 
graphitic regions. 

Existing theories of polishing, based on microcleavage 
B, are, as far as the soft direction is concerned, not con- 
sistent with this new experimental evidence. In the hard 
direction, this experimental evidence does not suggest 
cleavage along the (111) planes but does not fully rule 
out a different fracture mechanism. As yet no detailed 
atomic mechanism explaining the polishing anisotropy 
has been proposed. 

To investigate the asymmetry in the wear rate and the 
differences in post-polishing surface morphology we per- 
formed simulations in which a single nano-asperity, con- 
structed on the (110) surface, was deformed by a rigid 
tip. We first present the results of separate simulations 
in which the tip was incident from the soft and the hard 
polishing directions, which show pronounced differences 
in the extent of the induced deformation. Then these 
differences in asperity removal are related to the pro- 
cess of nano-grooving which, we believe, is responsible 
for the polishing anisotropy. Finally we present a simple 
schematic representation of the nano-grooving process 
which highlights the qualitative differences between the 
directions. This wear process is predominantly mechani- 
cal rather than chemical in agreement with the suggestion 
of Couto et al. [g. 

The implementation of a standard total-energy pseu- 
dopotential method JTo| on a massively parallel computer 
permits the quantum mechanical simulation (with an ac- 
curate description of bond breaking and formation) of 
systems containing several hundred atoms, sufficient for 
the modelling of a realistic nano-asperity and tip. 

The asperity was constructed by adding tetrahedrally 
bonded carbon atoms onto a clean, unreconstructed, 
(110) surface. All dangling back-bonds were then sat- 
urated with hydrogen and the system was allowed to re- 
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lax. The base plane of carbons of this relaxed system 
was then held rigid during the subsequent simulations. 
We chose to saturate the surface with hydrogen because 
of the large increase in polishing rate in hydrogen-rich 
environments and the known affinity of diamond surfaces 
for hydrogen jflj. We formed a rigid tip, with which to 
deform the asperity, from a fragment of a Pandey (2x1) 
reconstructed (111) surface (l^JT^]. The dangling C-C 
bonds of this fragment were saturated with hydrogen. 
The (111) surface is realistically modelled as a rigid tip 
as this surface is the natural cleavage plane of diamond 
and is known to be relatively hard in all directions. The 
Pandey reconstructed tip was chosen in preference to a 
hydrogen saturated (111) surface as, although polishing 
is performed at temperatures below the desorption tem- 
perature of hydrogen, the wear process will remove any 
hydrogen from the abrading diamond. For both direc- 
tions the rigid tip was advanced into the asperity in steps 
of 0.1 A. At each position we allowed full relaxation of the 
remaining atomic degrees of freedom. This adiabatic ap- 
proximation is justified as the speed of sound in diamond 
is approximately three orders of magnitude greater than 
the typical velocity of the abrading diamond. Figure |l| 
shows the initial atomic configurations of the supercells 

0- 

The local density approximation was used to describe 
exchange and correlation, and the electronic states were 
expanded at the T point of the Brillouin zone using a 
plane-wave basis set truncated at 44 Ry. The calcula- 
tion was converged with respect to plane wave cutoff and 
supercell height to an accuracy of O.leV per atom. The 
electrostatic interaction between the asperities in neigh- 
bouring supercells was monitored during the simulations 
and was found to be insignificant. An optimized non- 
local pseudopotential jL{| was used to describe scatter- 
ing from the Is core of the carbon atoms and a Coulomb 
potential was used for hydrogen. The non-local pseu- 
dopotential was applied in the Kleinman-Bylander form 

s 

The force felt by the tip as a function of its position 
(Figure [^) shows some obvious differences between the 
hard and soft directions. In the hard direction, the tip 
can be advanced much further before reaching a yield 
point, beyond which all rigidity is lost. Hence, since the 
forces on the tip are of comparable size for both direc- 
tions, the impulse delivered to the tip is larger in the 
hard direction. However, this is not sufficient to explain 
the anisotropy in the polishing rate. Another feature of 
the hard direction is the large spike in the forces which 
might deflect the incident tip. 

The reason that the normal force is larger in the soft 
direction is that the tip was at a slightly lower initial 
height above the surface during this simulation. The 
height and angle of the tip were chosen to make rea- 
sonable contact with the asperity for the simulations but 
this is not believed to be critical to the mechanism. The 



justification for this assertion is evidence from simula- 
tions in which the rigid tip was replaced by a term in 
the Hamiltonian representing fictitious mechanical forces 
p7[ . These forces had no component normal to the sur- 
face and the asperity was allowed to relax as the forces 
were steadily increased. The results of these simulations 
indicate that the process is predominantly mechanical 
and that the role of the normal force is not important 
to the mechanism in the soft direction although it does 
help to stabilize the asperity in the hard direction by pre- 
venting it from peeling off. In view of this evidence we 
expect no qualitative change in the mechanism provided 
that the tip makes reasonable contact with the asperity 
but does not cause direct subsurface damage. 

Figure || is of particular interest because it allows us 
to identify the particular atomic configurations which are 
responsible for changes in the asperity's degree of resis- 
tance to deformation. Points where the key processes 
occur have been labelled in the figure and the configu- 
rations corresponding to these points are shown in Fig- 
ures H & |]. Points relating to the soft direction will be 
discussed first. At position (A) both the normal and 
the retarding forces have been slightly reduced by the 
loss of a hydrogen from the asperity and the subsequent 
formation of a C-C bond between the asperity and the 
tip |l8). This rearrangement does not alter the tetra- 
hedral coordination of the atoms in the asperity which 
persists until the first C-C bond in the asperity is bro- 
ken (B). Beyond this point the asperity is compressed, 
mainly through bond angle deformation, until, between 
(C) and (D), there is a reconstruction in which the as- 
perity loses its remaining structural strength through the 
breaking of a six-membered ring. At this stage the ab- 
sence of structural rigidity permits the few remaining 
single C-C bonds, by which the asperity now adheres 
to the surface, to be broken with relative ease. At the 
termination point, the structure of the asperity is neither 
diamond-like nor graphitic and any polishing residue is 
likely to be amorphous |l9| . Throughout the process the 
deformation remains localized, affecting only the atoms 
within the asperity itself. This suggests that the surface 
left behind will be flat on the atomic scale. A quanti- 
tative estimate of the upper limit on the wear rate was 
made using numerical data from the simulation, under 
the assumption that this process is entirely responsible 
for the experimentally measured normal load jf7||. This 
estimate is significantly above the maximum experimen- 
tally observed wear rate indicating that the process can 
easily account for the magnitude of the observed wear. 

Turning now to the hard direction (Figure |J), posi- 
tion (E) is where a carbon in the tip, which had already 
bonded to a hydrogen atom from the asperity, breaks 
this bond, and rebonds to a carbon atom in the asperity, 
leaving the hydrogen free to bond to a different carbon in 
the tip. As the tip is advanced towards (F) more bonds 
are formed between the tip and the asperity, although 
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in supplying the large force on the tip the asperity has 
only deformed to a limited extent. Between (F) and (G) 
the force on the tip falls off dramatically and the region 
of the asperity in the immediate neighbourhood of the 
tip reconstructs. The effect of the reconstruction is to 
force the two carbons at the leading edge of the asperity 
beneath the tip. This causes significant sub-surface dam- 
age which may even have been slightly restricted by the 
rigid base plane. The structure after the transformation 
looks similar to the structure at the position (E) where 
the force began to accumulate. This suggests that, if 
the asperity were longer, and therefore more able to sup- 
port a large retarding force, the forcing of the carbon 
atoms into the sub-surface could be repeated. At (H) 
the remaining asperity carbons are in a highly distorted 
configuration, primarily due to large bond angle distor- 
tions, and the atoms at the back of the asperity yield first 
leading to a complete loss of rigidity. At point (I) the re- 
tarding force has fallen significantly and it is likely that, 
as in the soft direction, the remains of the asperity would 
form an amorphous residue. However, in contrast to the 
soft direction, the deformation in the hard direction is 
highly non-local. 

The process of removing an asperity may be related to 
that of forming a nano-groove by making two changes to 
the system. Firstly, a change in boundary conditions at 
the edge of the asperity is required to recreate an ideal 
surface. Secondly, the chemical environment beneath the 
tip must be altered owing to the removal of the preced- 
ing asperity. The effect of making these changes is very 
different for the two directions. In the soft direction the 
change in boundary conditions involves linking the (110) 
chains to join asperities in adjacent supercells. A change 
in the chemical environment beneath the tip will mainly 
affect the normal force which, as already discussed, is not 
believed to be important. The local nature of the defor- 
mation during asperity removal suggests that both the 
mechanism and the forces occurring during nano-groove 
formation will be similar to those observed in the simu- 
lations. In contrast, in the hard direction the forces oc- 
curring during nano-groove formation are expected to be 
significantly higher than those arising during the removal 
of an asperity, although the mechanism will be similar. 
The change in boundary conditions has the effect of ex- 
tending the rear of the asperity to infinity. Consequently 
it is able to supply much greater resistance to the on- 
coming tip. The chemical environment beneath the tip 
may be significantly altered as the sub-surface is likely to 
be damaged. As already mentioned, the process of forc- 
ing carbon into the sub-surface may be repeated until 
either the increasing normal force due to the compressed 
carbon deflects the tip or bond breaking occurs at an in- 
termediate distance from the tip leading to the removal 
of carbon. Figure |^ shows the suggested nano-grooving 
mechanisms. 

In summary, we propose the mechanical process of 



nano-grooving as the principal mechanism for wear of 
the (110) surface of diamond during polishing in the soft 
direction. The mechanism proceeds mainly via bond an- 
gle deformation, and stress is localized on a few diamond 
bonds at a time, allowing bond breaking. The role of the 
normal force does not appear to be important, the polish- 
ing residue is expected to be amorphous and the surface 
will be left atomically flat. In the hard direction nano- 
grooving may play a role in the wear of diamond but the 
process is significantly more difficult. Provided that a 
normal force is present to stabilize the surface the defor- 
mation is not local and the stress is distributed over a 
large number of diamond bonds. There is extensive sub- 
surface damage and the likely polishing residue would 
again be amorphous. We believe that these differences 
in the atomic mechanisms of nano-groove formation are 
responsible for the observed anisotropy in the polishing 
rate. 
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We define a bond to be formed when the relevant atoms 
are separated by less than the sum of their covalent radii 
i.e. 1.72 A for a C-C bond and 1.22 A for a C-H bond. 
Beyond this point, the absence of rigidity made the Born- 
Oppenheimer surface shallow and hence the validity of 
the adiabatic approximation became questionable and 
the calculation became expensive to converge. 
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FIG. 1. Relaxed initial atomic positions: (Left) Soft direc- 
tion simulation viewed along (110). This supercell has a vol- 
ume of 1291 A 3 and contains 93 C and 68 H atoms. (Right) 
Hard direction simulation viewed along (001). This super- 
cell has a volume of 1603 A 3 and contains 111 C and 80 H 
atoms. (Center) Top view of the (110) surface showing the 
two scanned directions: (110) parallel to the characteristic 
chains and (001) perpendicular to them. 
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FIG. 2. Forces on tip during the simulation. Open symbols 
refer to the soft direction, filled to the hard direction. Trian- 
gles refer to the retarding force, i.e. the force which opposes 
the advance of the tip, and circles refer to the normal force. 



FIG. 3. Atomic configurations, corresponding to the points 
shown on the force graph (Figure ^), relating to the simula- 
tion performed in the soft polishing direction. The diagrams 
correspond to tip positions (A) = 0.6 A, (B) = 1.2 A, (C) = 
1.8 A and (D) = 1.9 A. 




FIG. 4. Atomic configurations, corresponding to points 
shown on the force graph (Figure ^), relating to the simula- 
tion performed in the hard polishing direction. The diagrams 
correspond to tip positions (E) = 1.7 A, (F) = 2.4 A, (G) = 
2.7 A, (H) = 4.4 A and (I) = 4.9 A. 
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FIG. 5. Carbon skeleton schematic representation of sug- 
gested nano-grooving mechanism. The top figure is for the 
hard direction and shows the deformation to be non-local. 
The lower figure is for the soft direction and shows the mech- 
anism to be local and repeatable. 
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